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Abstract 

We investigate the three-dimensional SU(3) gauge theory at finite 
temperature in the framework of dimensional reduction. The large scale 
properties of this theory are expected to be conceptually more compli- 
cated than in four dimensions. The dimensionally reduced action is com- 
puted in closed analytical form. The resulting effective two-dimensional 
theory is studied numerically both in the electric and magnetic sector. 
We find that dimensional reduction works excellently down to tempera- 
tures of 1.5 times the deconfinement phase transition temperature and 
even on rather short length scales. We obtain strong evidence that for 
QCD 3 , even at high temperature the colour averaged potential is rep- 
resented by the exchange of a single state, at variance with the usual 
Debye screening picture involving a pair of electric gluons. 
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1 Introduction 



Dimensional reduction is a powerful technique to study field theories at high temper- 
ature. In the Euclidean formulation such theories are defined in a volume with one 
compact dimension of extent 1/T. As the temperature T becomes large, one may 
expect that the non-static modes in the temperature direction can be neglected, and 
that one is left with a theory in one less dimension [p], This naive reduction is, 
however, only exact in the classical limit. In general, one has to integrate over the 
nonstatic modes to obtain the effective action |3|, f|, || . As was shown in Refs ||. [| 
the effective action for the long distance phenomena will, however only contain a 
limited number of local terms at sufficiently high temperature, because the integral 
over the nonstatic modes does not contain infrared divergencies. Furthermore the 
coefficients can be determined from a perturbative expansion of the integral over 
the non-static modes. This effective action is then expected to describe correctly 
the infrared behaviour of the full theory. 

This method was first successfully applied to lattice gauge theories in Refs || 
H 0]. In particular the screening mass as defined from the correlation between 
Polyakov loops, which describes the colour averaged potential at finite temperature, 
was calculated. The renormalized coupling constant, used in the reduction, was 
defined from a perturbative calculation of this potential at small distances. In the 
application to pure lattice gauge theories it was possible to compare the screening 
length in the reduced theory with that calculated in the full theory. It was found 
that the dimensional reduction gave good results in the gluon plasma phase down 
to T/T c w 2, where T c is the transition temperature between the confined and the 
plasma phase. 

The real challenge is of course the calculation in full QCD. Here the fermions 
can be integrated out explicitely, because they have no static modes, and they only 
change the coefficients in the effective action. A calculation of the screening lengths 
was performed in ||. A summary can be found in II. 

Later the formalism was adapted to and successfully used to investigate the 
electroweak phase transition in [[Kj. It was also employed to determine the chro- 
moelectric screeening mass in QCD over a large range of values of the temperature 
TT|, O, |T3|. In these articles another operator was used to define a screening mass. 



Also another renormalization scheme was used, the MS scheme, leading to a renor- 
malized squared coupling around three times larger in the similar temperature range. 
One may thus expect higher order corrections to be more important. Further evi- 
dence for dimensional reduction in the region of a few times T c have been found in 
I3J, [0 and m. 



It is, of course, very interesting if several quantities can be calculated in the 
same scheme. In particular quantities depending on the chromomagnetic sector, 



2 



which have infrared singularities in perturbation theory that cannot be determined 
so far by resummation, need a non-perturbative approach. One gauge independent 
operator, which has been studied is the spatial string tension. However, it is still 
an open question what the corresponding observable is in the reduced model. In 



Ref. JlTj it was shown that this string tension increases with the temperature in the 



deconfined phase. Later it was compared with the three-dimensional naively reduced 
theory, and a good agreement was found |18|], using, however, the bare coupling 
constant, which even on moderately sized lattice is different from the renormalized 
one employed in Ref. || - ||. 

In the present article we want to show more precisely how well dimensional reduc- 
tion works, both for correlations between Polyakov loops and for the spatial string 
tension, in the case of three dimensional SU(3) gauge theory at high temperature. 
The effective model in this case is a two-dimensional Higgs model, with the Higgs 
field in the adjoint representation. The interest of this calculation is that the three- 
dimensional theory has many properties in common with the four-dimensional one, 
but the numerical simulations become much simpler, and one can obtain much more 
precise information on the quality of dimensional reduction. However, the infrared 
behaviour of perturbation theory is worse, so that even the leading perturbative 
definition of the Debye mass is not well defined. Three dimensional gauge theories 



at finite temperature have been considered earlier by D'Hoker [19 



In this article the three-dimensional model is investigated through numerical 
lattice simulations at finite temperature. Dimensional reduction is performed by 
perturbation theory in the static time averaged Landau gauge, using the lattice reg- 
ularization, when necessary. The effective two dimensional model is also investigated 
by lattice simulations, which automatically include the non perturbative features of 
the model. A detailed comparison of correlations between Polyakov loops and of 
the spatial string tension is presented. Because the spatial dimension is two, we 
expect that the infrared properties of the theory are even more complicated than 
in 3 + 1 dimensions. In particular, the Debye screening picture does not necessarily 
translate to QCD 3 . 

Our paper is organized as follows. In the next section we introduce the lattice 
implementation of 3-dimensional SU(3) QCD at high temperature and derive the 
associated dimensionally reduced model. Special care is taken in defining the self 
couplings of the adjoint scalar of the 2D action in the limits appropriate to the 
large temperature behaviour of QCD 3 in the continuum (scaling limit). Section 3 
is devoted to the corresponding numerical simulations in 2 and 3 dimensions, and 
to the measurements of Polyakov loop correlations and spatial Wilson loops. The 
results are presented and discussed in Section 4. We show that the Polyakov loop 
correlations measured in (2+l)D or through their 2D reduction agree at large dis- 
tances, and remain close to each other down to quite small distances, the more so the 
temperature is high. This finding gives quantitative support to the statement that, 
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for static observables, the reduced action correctly describes the (2+l)D properties 
in the limit where the momenta are small compared to temperature. The decay of 
these correlations with spatial distance is not compatible with the Debye screening 
expected from the exchange of two electric gluons, but strongly favours the existence 
of an isolated excitation whose mass is measured. We also show that, although the 
Wilson loops are non-static operators, the (spatial) string tension in (2+l)D is close 
to that measured in 2D, which itself differs only weakly from the known pure gauge 
2D string tension. We summarize and conclude in a last section. 



2 Lattice Action in 3D and its Reduction to 2D 
2.1 Lattice Definitions 

In this subsection we give the definitions pertinent to the simulation of SU(3) gauge 
theory in 2 + 1 dimensions, at finite temperature. 

The gauge field action is taken to be the Wilson action on a cubic lattice A 3 , 
namely 

2 

S^ v ( y U)=p 3 J2 Yl ^-^^U^^Uix + a^^Uix + a^fi^Uix-^)- 1 ^. (1) 

The axis ji = is along the Euclidean time or temperature direction, and a is the 
lattice spacing. An SU(3) gauge group element U (x; //) is related to the gauge fields 
Afa),de [1,8], by 

U(x;n)=exp[iA li (x)], (2) 
A,{x) = j^Al{x)\ d , tv\ d \* = U de . (3) 

d=l 

The A's denote the 577(3) algebra generators. The gauge fields are dimensionless. 
The lattice coupling f3 3 is related to the conventional bare gauge coupling g 3 of 
perturbative QCD by 

ft = —2- (4) 

The theory is super-renormalizable, g%, which has dimension one in energy, is pro- 
portional to the renormalized coupling constant squared, and can thus be used to 
set the physical scale. The 3D lattice action describing the system at a finite tem- 
perature T is obtained by taking for A 3 a cubic lattice L 2 S x L , and 
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The observables of interest in this work are the spatial string tension, as extracted 
from the expectation values of spatial Wilson loops formed by links i joining the 
sites x = x + xo of a rectangle R = [R\, R%] at xq fixed, 

W(R 1 ,R 2 ) = hrl[U(x-?), (6) 

and the expectation values of Polyakov loops L(x) and their connected correlations 
P(x) 

L{x) = -tr TJ U{x + rr aO,0). (7) 
■J , 

x =l 

P{x) = (L{c)L\c + x))-\{L{c))\ 2 - (8) 



We want to study the behaviour at large T of the continuum limit of the theory 
defined by S^- Before proceeding, let us state more precisely what these continuum 
and large T limits mean in terms of the lattice parameters L , 0% and L s . Given g 3 
and T we have from 



6 AT 1 

— , and L = — . 
agl a 1 



P3 = 7-3, and L = —. (9) 



It follows that as a — > 0, scaling (constant physics) corresponds to 



]h_ = T 
6L gl 



00, /3 3 — > 00, t = -^-j— = — = constant. (10) 



The dimensionless quantity r thus measures temperature in units of the scale gf , 
high temperature means large r values. From numerical simulations we obtain that 
T c j g\ ~ 0.61 [20 1 . Also, on the lattice, L s must be kept much larger than the largest 



spatial correlation length in lattice units. Assuming that this length is of order 
l/( a ^3^^ 1_ ) f° r some oc > 0, this together with ([5]) implies 




(11) 

Note that this ratio increases with temperature. As we shall see, these consider- 
ations, important for the design a meaningful numerical investigation of the high 
temperature phase of QCD 3 , are also useful to discuss the adequacy of its 2D re- 
duction to describe at least part of its large distance physics. 
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2.2 Dimensional Reduction: the Effective 2D Action 



The dimensional reduction means that one should derive an effective action for the 
static modes, obtained by keeping the static part of the action and adding terms 
coming from a perturbative integration over the non-static modes. It is therefore 
convenient to choose a static gauge, where this separation can be made. 

The operators of interest to us are gauge invariant and may be computed in 
any gauge. In particular, a static gauge can be realized on the lattice with periodic 
boundary conditions by choosing A independent of the imaginary time coordinate 
x : 

A (x ,x)=A (x). (12) 

In such a gauge the Polyakov loops are static operators, unlike spatial Wilson loops, 
and read 

L(x) = -tr exp[z L A (x)]. (13) 
3 

Full gauge fixing, which is necessary for forthcoming perturbative calculations 



211 , is achieved by adding the Landau constraint 

2 

Y} A ^ - A ^ x - ^ = °- ( 14 ) 

xq i=l 

Together with ([12]), this condition constitutes the so-called static, time-averaged 
Landau gauge (STALG, see [^, |23|j). 



At high temperature and in the thermodynamical limit, the well known Z3 sym- 
metry of the action is broken (deconfined phase) and the average Polyakov loop 
(L(y)) is not zero. Its phase may be either one of the three 3 rd roots of one, the 
corresponding threefold degeneracy can be lifted up by rotating this phase away, so 
that the U's are always connected to the unit matrix. 

In what follows, then, the gauge manifold will be parametrized by gauge fields 



A Ax) subjected to the constraints ( |12] , [L4|), and perturbation theory, whenever it 



applies, is generated by expanding the 3D action into small fields. 

The effective 2-dimensional action is obtained according to standard techniques 
^ consists in the classically reduced 3D action, with all fields restricted to 
their static components (no dependence in the th coordinate), supplemented by the 
interactions generated by integrating over the non-static degrees of freedom which 
are left over after gauge fixing, namely the non-static part A^ s (xq,x), i=l,2, of the 
gauge fields, constrained by (fH]). The resulting action depends on the 2-dimensional 
SU(3) gauge fields Af aUc (x) defined by 

A t (x) = Af aUc (x) + A?(x ,g), (15) 
J2Ar(x ,x)=0, (16) 

XO 
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and on Aq(x), a scalar in 2 dimensions (Higgs field in the adjoint representation). 
In the following, the superscript static on the 2D gauge fields will be omitted. So 
defined, this action is of course non polynomial and non local in A . 

At large T however, one may truncate its expansion in powers of A(x), adapting 
the method developed in details in || for the 4D^3D reduction. In the present 
case, power counting and Becchi-Rouet-Stora symmetry of the gauge fixed theory 
leads to retain only the 1-loop contribution of the non-static modes to the quadratic 
and quartic A Q self-couplings at external momenta p = 0. The contributions of 
monomials of higher degree, or non local in A Q are suppressed by powers of (g 2 /T) 2 
and (p/T) 2 , up to log's. More precise statements will be made at the end of this 
section. 



We write down the explicit result S 2 ^ of this reduction, in terms of the param- 
eters L and /3 3 of the initial 3D lattice action (L s is assumed to satisfy (O) so that 
finite size effects are negligible). 

For later convenience, we rescale the scalar field by setting 



4,(aO = ^^JL, (17) 

and express S 2 jj as a function of <fi and of the 2-dimensional U's. It contains three 
pieces, a pure gauge part S^, the gauge covariant kinetic term Su,4> for 0, and the 
self-interaction induced by integration over the non static fields: 

S 2 ff (U,<f ) ) = S 2 v + S u ^ + S ( t > . (18) 

Here S^y is the Wilson action for D = 2 obtained by repeating L times the purely 
spatial part of S^r. Hence 

= fa ^(l--^tiU{x]l)U{x + a%2)U{x + a%l)- 1 U{x]2)- 1 \ (19) 

zeA 2 

& = (3 3 L . (20) 

The lattice kinetic term for the (adjoint) field <fi is obtained by expanding the S^r 
action to second order in A , and reads 

2 2 

s u,<» = EE tr ( A( w (f) ) > ( 21 ) 

X 8=1 

Di(U)(f)(x) = ?7(f;z)0(f+aT)?7(f;z) _1 -0(f). (22) 
Finally the "Higgs potential" is written 

S <t> = ^h 2 ti<p{x) 2 + h i (ti4>(x) 2 ) 2 . (23) 
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In the scaling limit ( |TOD explained in subsection 2.1, we drop those terms of the 
small a expansion of h 2 /a 2 and h^/a 2 which vanish as a —>■ and find 

h > = vk%> (24 > 

9 / , , \ . 5 

h 2 = 



7T L fa 



■(cq log L + a ) ; c = 1, ci = -log2-l. (25) 



We call these expressions the scaling parts of h 2 and /i 4 . Let us sketch how these 
results have been derived and comment on them. 

The coupling requires the 1-loop nonstatic contribution to the Green's func- 
tion Tn)(p = 0) with four external A fields at zero momentum. At T = 0, the 
Feynman integrals involved in the continuum all are of the type J d 3 k/k 4 and thus 
UV convergent but IR divergent. At T ^ 0, the restriction to non static modes 
(k 7^ 0) makes them IR finite. For these reasons, can be, and was, calculated 
directly in the continuum. Our result (]24]) coincides with that obtained by adapting 
to D = 3 the calculation of Nadkarni || for D = 4. 

The case for h 2 is slightly more tricky. It is given by 

h 2 = -^n-(O), (26) 

where I1qq(0) is the 1-loop contribution of the non static modes to the vacuum 
polarization at zero momentum. Two Feynman graphs contribute, with one and two 
Ai internal propagators, respectively of the form g 2 J d 3 k/k 2 and g% J d 3 k /cq/(/c 2 ) 2 , 
both linearly divergent with the cut-off A. However, at T = 0, gauge invariance 
guarantees n o(0) = 0, so that the two graphs cancel each other exactly there. At 
T^0, the cancellation is only partial, still leaving a logarithmic UV divergence. We 
thus computed I1qq(0) on the lattice, which provides a cut-off of order A ~ 1/a — 
T L at high temperature according to Eq.@, hence the log(L ) term in (p5|) . 

Applying the Feynman rules of f24]| , the finite lattice expression of I1qq(0) is 
found to be 

s k ^0 £ ( k 2 \ K 
2 

^ = 2sin(V2), k 2 = Y,^\ % = cos(V2)- (28) 

o 

The L s — > oo limit is obtained by the replacement (27r/L s ) 2 ^^ — > j^ n d 2 k. After 
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some manipulations, one arrives at 



nSJ(O) = ±^(-l(l-2t)S(t)-t(l-t)*?&), i = S m 2 (^), (29) 



1 [* dkldkl 
{t) = 4^ J_ n 2(1 + t) - cos(fci) - cos{k 2 ) ' (30) 

It is clear on the latter equation that S(t) has a logarithmic singularity at t = 0. 
The result ( p5|) follows from (^) after the asymptotic behaviour at large L and 
infinite L s of I1qq(0) is inserted. Details on its derivation from (|29|) are given in the 
appendix. Here, we just mention a useful trick used to evaluate I1qq(0) , noting that 
Eq.(^) can be rewritten as 

BSC) = -^£>M- * = T?. CD 

n =l 

= sina; cosa; 5( sin 2 x). (32) 

At large L , En =i £ G ( X ) is of order A)A f;~" /io rfG(x) = -2L /irG(n/ L ): The 
linear divergence in L is cancelled by the sin(7r/L ) factor in G while the logL 
originates from the logarithmic behaviour of S^sin 2 ^) at x ~ 7t/Lq. 



2.3 About the Reduced Action: Summary and Remarks 

We first summarize what has been done above. The QCD3 lattice action is given 
by Eq. ([!]). The parameters of physical relevance are the lattice coupling @ 3 and L . 
The observables we consider are the Polyakov loops (0) and their correlation (|3|), 
and the spatial Wilson loops (^). 

Dimensional reduction leads to the 2D lattice effective action described by Eqs. 
([T8| - |25D . It is a gauge invariant model for a scalar in the SU(3) adjoint representa- 
tion, whose "Higgs potential" is truncated at order 4 in at high temperature and 
low momenta. The lattice gauge coupling and scalar self couplings 02, ^2 and /14 are 



fixed by the 3D parameters Eqs. (^, [25|, ^). The Higgs field is normalized in such 
a way its kinetic term in the classical limit is tr (<9j0) 2 . 

Let us now comment upon what we have obtained. First, as announced before, 
it is consistent to truncate the effective action to order 4 , with I12 and couplings 
computed at 1-loop order. To be more precise, power counting and inspection show 
that h n /a 2 , with h n the coupling associated with any monomial of degree n in 0, 
has the following order of magnitude: 



n 



g\T + g\ T logT (l + 0[(p/T) 2 , gl/T]) , (33) 
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n = 4 
n > 6 even 
n > 3 odd 



<? 3 4 (1 + 0(^3 2 /T)), (34) 
(T 2 (sg/T) B/a ) , (35) 
(T 2 (rf/r) B/a (P/T)) • (36) 



The corrections for n = 2 and 4 come from 2-loop order and non-local effective 
interactions. It is thus consistent to neglect the n = 6 and higher power monomials, 
as well as all couplings with n odd, and to keep h 2 as given by Eq. (|2~5f) . 

We next note that unlike h± which is positive, insuring boundedness of the parti- 
tion function at large fields, the quadratic coupling h 2 is negative. The "potential" 
thus have a shape typical for gauge symmetry breaking by Higgs mechanism. Of 
course, this symmetry breaking is not expected to occur in 2 dimensions. In fact, 



the form (p25|) of this coupling can be seen as a counterterm for the logarithmic UV 
divergence of the 2D model, here prescribed in advance by the UV regularization of 
the 3D model, while the actual IR behaviour of the theory is highly non perturbative. 

In the static gauge chosen, the 3D Polyakov loop operator is given in 2 and 3 
dimensions by the same function ([13]) of the Aq field. According to the normalization 
(JTTP of the field and to Eqs. Ql0|, fZ0|), it reads 

1 



x 



tr exp(z </>(£) /y/r), (37) 



ft ft T 



6L 6L 2 , gl ' 



(3* 



Along a line r = const, (constant physics, see subsection 2.1), the operator used to 
probe dimensional reduction remains unchanged, and it is thus consistent to compare 
the Polyakov loop correlations measured in 3 and 2 dimensions at the same fixed 
r value. The continuum limit is L — > oo, and constant (large distance) physics at 
large L can be checked within either one of the two models. 

Let us then consider the effective action S^t in terms of the dimensionless tem- 
perature r, Eq. d38|) , and of the lattice parameter L . After setting a 2 g 2 = 6/ft, 
and using Eqs. (|20| , [25| , [24| , p8] ), the three dimensionless couplings (a^) 2 , h 2 and h± 
read 

M 2 = -L (39) 

TLln 



3 c log L + a 
2tt tL\ 



(40) 



111 

647T tLq t 

We see that, except for the logL in h 2 , these three couplings all scale as 1/Lq with 
Lq, while as a function of r, h^ is r times smaller than the two other ones. 
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Let us finally write down the effective Lagrangian C e ff, as obtained from the 



small a expansion of the effective action S 2 ^. 



In S^ff, we make the substitution 



Ai — > ag 2 Ai, replace J2s by a ~ 2 / ^ 2;r ' use and g\ as parameters instead of the 
lattice /?3 and L , and take the limit a — > (but in logL = — log aT in the <fi 2 
term). We obtain 



C 



eff 



Di 



F 



C 



CT 



j£^5 + *[^] 2 + J;(§)W + z: 

c=l 

di(p + ig 2 [Ai,(/)], 
diAj - djAi + ig 2 [A i ,A j ], 
3<? 2 2 



CT, 



(42) 



2tt L 



-log(aT) + 5/2 log 2 



tr 



(43) 



This is a 2D, SU(3) gauge invariant Lagrangian for an adjoint scalar 0, but it 
is far from being the most general one. The gauge coupling g 2 , with its canonical 
dimension one in energy, sets the scale. The non kinetic quadratic term is the 
counterterm Cct, suited to a lattice UV regularization with spacing a. The quartic 
self interaction A4 = g\ x (g 2 /T) 2 , generically a free parameter, here goes to zero 
as T — > 00 in units of g 2 , and presumably plays a marginal role. We also note 
the absence in the model of odd powers <p 2k+1 , allowed by gauge symmetry for 
k > 0. The corresponding Z 2 symmetry might however be spontaneously broken in 
some subspace of the unrestricted parameter space g 2 ,h 2 ,h^. The similar problem 
in 4 — > 3 QCD reduction has been studied in Ref. (T^j. For the case of SU(2) 
in (3+l)D, Z 2 is the center of the gauge group so that Z 2 breaking is also gauge 



symmetry breaking, a subject previously discussed in Refs. | 25| , [ 26[| and [11 



3 The Numerical Simulations in (2+1) and 2 Di- 
mensions 



We simulated the reduced theory using the same method as in reference [0. For 
updating the gauge fields U we employed a mixed algorithm, where new trial SU(3) 
matrices were generated by the heath-bath algorithm and are accepted according 
to a Metropolis condition on the hopping term. We used the standard Cabibbo- 
Marinari pseudo heath-bath algorithm based on updating the SU(2) subgroups of a 
SU(3) matrix [p^] . The new SU(2) sub-matrices were generated using the Kennedy- 
Pendleton algorithm [p8|| . After each subgroup update the resulting matrix was 
subjected to the Metropolis question on the hopping term. This approach results in 
an acceptance rate close to 95%. 
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Alternatively we tried the multi-hit Metropolis algorithm. An update of a single 
link variable was attempted several times (typically eight) before attempting to up- 
date next variable. The trial matrices were generated according to : U — > exp(iA)U 
where A was an hermitian matrix generated from the distribution exp(— e trA 2 ). 
The parameter e was tuned to obtain the acceptance ratio around 50%. Because 
generation of trial matrices is expensive in CPU time (matrix exponentiation) we 
took the values from a table that was regenerated every five sweeps. The autocor- 
relation times for this algorithm were comparable to the heath-bath but the CPU 
time required to perform a metropolis update was approximately two times bigger. 

The scalar fields were updated using the multi-hit metropolis. At every point 
we tried typically eight updates. New scalar fields were obtained from A —>■ A + 
5 A where 5 A was generated according to exp(— e tiA 2 ). For each set of coupling 
constants, e was tuned to obtain the acceptance ratio around 50%. One sweep over 
the lattice consisted of one update of gauge fields followed by one update of scalar 
fields. After each sweep we measured "control" variables : expectation value of 
plaquette, tr A 2 , tr A 3 and tr A 4 . 

All runs except for scaling test (see the next section) were performed on 32 x 32 
lattices with the parameter L set equal to 4. The dynamics of the gauge and Higgs 
sectors turned out to be quite different. The integrated autocorrelation time of 
the plaquette operator was typically non-measurable (< 1) while for tr A 2 it varied 
between 100 and 400 sweeps when /3 3 was varied from 21.0 to 346.0. For each value 
of P we have collected typically 10 6 sweeps. 

In three dimensions we used the data previously gathered on the Quadrics APE 
supercomputer p0| on 32 x 32 x 4 lattices. 

The main quantities measured were the average (L(x)) of the Polyakov loops, 
their on-axis 2-body correlations -P(r), and the average (W(Ri, Rq)) of spatial Wil- 
son loops. In three dimensions they are defined by (y), (0) and (§). The only 
difference in two dimension is that the Polyakov loop is given by ([13]). When neces- 
sary, an index D = 2 + 1 or 2 specifies the dimension considered. Both 2D and 3D 
data were analyzed in the same way as to minimize the systematic errors. 

In 2D simulations Wilson loops were measured every ten sweeps and blocked by 
50. Integral autocorrelation time for the measurements of 14 x 14 Wilson loops was, 
in the worst case, of the order of one and was negligible in the majority of cases. 
The measurement routine is quite expensive and consumed almost half of the CPU 
time. 

In order to extract the string tension, local potentials were first extracted from 
the ratios : 
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By definition, the potential is 

V(R) = lim V(R,&). (45) 

R'— >oo 

For each given R, V(R, R') was found to decay exponentially to a constant in R' . In 
practice, the constant regime is reached within errors above R' = 3 and the potential 
was fitted in the range R' £ [4, 12]. The string tension a was then obtained from 
the ansatz 

V (R) = V + crR. (46) 

A fit to the above linear potential is sufficient to obtain a stable value of a within 
errors when it is performed in the range from R m i n = 2 to R ma x- The R ma x values, 
depending upon the value of /3 3 , varied from 10 to 14. The \ 2 value in the 3D case 
was of the order of one, and about ten times less in the 2D case, likely due to large 
correlation between points. 

The errors on a were calculated by the triple application of the jack-knife method. 
In order to calculate the jack-knife estimate of the error on a from the fit to the 
formula (f46| ) we used 16 jack-knife copies of the potential V(r) with errors. To 
estimate the errors of every copy we applied the jack-knife algorithm with the same 
block size again to each of the corresponding 16 data sets, thus generating 15 copies 
of every set. However, in order to calculate the V(r) we needed the ratios V(R, R 1 ) 
again with errors. Calculation of these errors was done by a third and last jack-knife 
step applied to each of 16 x 15 data sets obtained in the previous step. 

Polyakov loops correlators were also measured every ten sweeps. While, due 
to the large auto-correlation time, this was in a sense an "oversampling" , it did 
not take much of the CPU time. The measurements were blocked by 50 and then 
written out. Afterwards the data were again blocked into 25 blocks. These 25 
blocks were used to calculate the connected correlators with errors estimated by 
jack-knife algorithm. We also tried different numbers of blocks and did not observe 
any significant change in the results. The double jack-knife algorithm, similar to 
the one described for the Spatial Wilson loops, was used to estimate the error for 
fits of the connected correlator to the formulae (HO, 150) used to extract a screening 



length (see next section). The fits were done in the interval r £ [r m i n = 4, 15]. We 
have checked that varying r min from 3 to 6 does not change the results. 



4 Results and Discussion 

We now present our results, discussing successively dimensional reduction (how well 
does the effective 2D model describe quantitatively the large distance physics of 
high temperature QCD3), scaling (is the parameter region explored close to the 
continuum limit), the screening length in the Polyakov loop channel and the spatial 
string tension. 
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4.1 Dimensional Reduction 



Let us compare the Polyakov loop correlations P2+1O") extracted from Ref. |20 
and P2O") measured in our simulation. A conservative statement on dimensional 
reduction is to say that the lowest physical state coupled to L(x) is the same in 
both cases, so that the two functions must have the same shape in r at large r. A 
stronger statement is that, to the extent that the weight associated with S^j- is a 
good approximation to the integral over the non-static fields of the (2 + 1)D weight, 
in the small p/T regime at least, the averages with respect to the two weights of a 
static operator such as (L(x) L(y)) are equal at large \x — y\. We will show that this 
latter situation is indeed approached at large enough temperature. 

We recall that temperature is given in units of the gauge coupling by 



gl 6L C 



Using L = 4 fixed, we investigate the temperature dependence of the correlations 
by varying (3$ in the range [21, 173]. High temperature means (3% sufficiently larger 
than the transition point in (2 + 1)D, which was found to occur at (3 C = 14.73 for 
L = 4 E(|. 



Fig.l compares the correlations P2+i( r ) and P^r) f° r ^3= 29 and 84. The curves 
are fits of the i^O") function, to be explained below. We see that the two correlations 
are extremely close one to the other: not only they have the same shape at large 
r, which is the primary prediction of dimensional reduction, but also nearly the 
same normalization. As announced, this finding favours dimensional reduction in 
the strong sense. In addition, we notice that this agreement between the two ways 
of computing the Polyakov loop correlation extends down to fairly small values of 
r. Since the normalization is set by \(L(x))\ 2 , it shows that, although it is local, the 
Polyakov operator is not very sensitive to the short wave length terms omitted in 
the effective action. That it is the case is born out by the data of Table 1, where our 
results for (L 2 (x)) and (L 2 +i(x)) are presented. Relatively to their distance to one 
(their common value at (3% = 00), their difference, already ~ 6% at /?3=21, decreases 
to ~ 0.6% at f3 3 =173. 

In Fig. 2, we illustrate more in details, for the (3$ values not reported on in 
Fig. 1, how well the 2D and (2 + 1)D correlations compare. Their ratio P2+1/P2 
is plotted against the distance r in lattice units. These data definitely support 
the statement that P2+1/P2 remains quite flat and close to one in the whole r and 
T range. Recalling that the lowest T value is only ~ 1.5 times T c , and that the 
correlation functions decrease in r by about three orders of magnitude, we conclude 
that the effective local 2D action reproduces the (2 + l)D Polyakov loop correlations 
with a remarkable accuracy, soon above the transition, and down to distances even 
shorter than 1/T. 
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Figure 1: Polyakov Loop Correlations Pd(j), D = 2 + 1 (circles) and 2 (squares), 
for /? 3 equal to 29 (left) and 84 (right). Distance r is given in lattice units. In both 
cases L s = 32 and L = 4. The 2D data were produced using h 2 from Eq. (p5J). 
The continuous and the dashed lines result respectively from fitting formulae 
and © to the 2D data. 



4.2 Scaling 

Given the temperature, the continuum limit is approached by taking 



T 



9l 



fixed, 



aT 



large. 



(48) 



Scaling is thus the statement that, at fixed r, the physics do not change with Lq, if 
Lq is large enough. We checked that by comparing the Polyakov loop correlations in 
2D for two sets of lattice parameters, namely [03, Lq, L s ] = [29, 4, 32] and [58, 8, 64]. 
Given T, doubling Lq means dividing the lattice spacing a by two, so that the 
physical spatial size aL s of the lattice is unchanged. The resulting correlations 
are presented in Fig. ^ showing a very similar shape as a function of the physical 
distance. It is again instructive to consider their ratio, found to be quite flat at 
all distances (Fig. |J): Within errors, there is no sizable deviation from scaling. 
Note that this constant ratio is not one. It is given by that of the corresponding 
values of \ (L(x))\ 2 , which is not a physical quantity (remember that the 2D effective 
Lagrangian contains an explicit (logarithmic) dependence on a via its counterterm 
Eq. 
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ft 


\(L2&)\ 


\(L2+i(x))\ 


21.0 


0.56002 (62) 


0.53467 (16) 


29.0 


0.67007 (25) 


0.66120 (13) 


42.0 


0.76397 (20) 


0.76130 (12) 


84.0 


0.87392 (11) 


0.87435 (6) 


173.0 


0.93494 (13) 


0.93530 (11) 



Table 1: The average Polyakov loop as a function of (3$ in (2+l)D and 2D. 

Hence scaling is verified in the range of interest. It justifies keeping Lq = 4, 
which is less expensive in computer time and allowed us to use the existing data of 
in (2 + 1)D. 

4.3 Screening Lengths 

The fast decay with r of the correlation functions is consistent with the existence of 
a finite spatial correlation length, £5, associated with the quantum numbers of the 
Polyakov loop operator, a colour singlet scalar. A mass can be defined as usual by 
Ms = £g X . In the similar situation with one more dimension, perturbation theory 
is invoked to state that Ms is 2Me, twice the so called electric screening mass, 
because the lowest "state" coupled to the loop is a two electric gluon state. If this 
is so, the P{r) correlation is expected to be proportional at large r to the square 
of the correlation function for a one particle state of mass Me- In the present 
case however, the infrared sector is more complicated and perturbation theory is 
doubtful. In particular we let open the possibility that there exists in the PL channel 
an independent screening mass Ms, associated with a true thermal excitation of the 
2D system. In the latter case, we parametrize the data according to 

^ = c (dw^' + K^-oK ^') ' (49) 

while in the former case one rather expects 

p( 2m E)( \ ~ r > ( \ p -m E r j \ -m E (L s -r)\ /rr,\ 

{T) ~ ° \[m E r}V* e + [m E (L s - r)K ) " ^ 

The m symbols denote masses in lattice units, i.e. m = aM. These parametrizations 
respect the lattice symmetry r — > L s — r. Since we have no direct access to tue, 
the two above expressions differ in shape through the prefactors, r -1 / 2 and r _1 
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Figure 2: Ratios P2+i(?")/-P2( r ) °f Polyakov Loop correlations as functions of the 
lattice distance r for /? 3 = 21,42, 173. 



respectively. We checked that for r > 1, lattice artefacts are negligible in the mass 
range considered. This was achieved by comparing Eq. Q4"9p to the lattice propagator 

PLatt(m s ,r) = — ^2 cos(pir)P Lati (m s ,p), (51) 

S Pl,P2 

P L att(™s,P) = 9 + 4sinh(m|/4). (52) 

in the notations of Eq. fl28|). We carefully analyzed our numerical data, and we 
find that the ansatz is by far the best one, giving a good fit of the data down 
to small r values. For illustration, the continuous curves of Fig. 1 are fits of the 
form (|49|) to P 2 (r),r > r min , with r min = 4, and changing r min from 3 to 6 does not 
change the answer for ms within errors. The only deviation which we observe is at 
(very) short distances, where some more massive contribution may be present, if not 
an artefact of the lattice UV regularization. It is anyway too small to be reliably 
analyzed. On the contrary, fits to the same data of the form Eq. flBTf ) lead in Fig. 1 
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Figure 3: Comparison of the Polyakov Loops Correlations in the 2D model for the 
two sets of parameters L , L s ] = [29, 4, 32] and [58, 8, 64] i.e. for constant r (Pp. 
The values of r for L = 8 (circles) are scaled down by a factor two in order to 
maintain the same physical scale. 



to the dashed curves, which are clearly not acceptable. The results of systematic fits 
of expression (49) to all available correlations P2O") an d -P(2+i)( r ) ar e presented in 
Fig.5. The corresponding dimensionless quantities Ms/(g?,\/T) = ms^/L /3^/6 are 
plotted versus g\jT = QL /j3s- Also shown for comparison are the masses resulting 
from a few simulations made with the so-called naive effective action (h 2 — /14 = 0). 
The latter are far away from the former, showing the important effect of taking into 
account the non static degrees of freedom. We also made runs where, instead of using 
the so-called scaling form (p5|) of h 2 , we took the value obtained numerically from 



Eqs. ([2g, ^Hp. The results are similar, especially at high temperature, but anyway 
distinguishable within our statistical accuracy. This shows a great sensitivity of the 
static properties to the quadratic counterterm (|43|) , an interesting feature per se. 
The scaling properties observed in subsection 4.2 of course reflect themselves 
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Figure 4: The ratio between the two Polyakov loop correlation functions presented 
in figure [3|. 

in the mass values. For the two cases compared there, we find M,sa=0.331(6) and 
0.169(3): the lattice spacing is reduced exactly by two within errors. 

Of course the fact that dimensional reduction works well for masses directly 
follows from its success for correlations. We preferred to illustrate it first on the 
latters, as done in Fig.2, that is independently of any interpretation of the nature 
of the observed screening lengths. What we now get in addition is that the main 
signal observed in the Polyakov loop channel favours the existence of a true colour 
singlet excitation of the high temperature gluon system, rather than it reflects the 
existence of an electric screening length (Debye screening) in the one longitudinal 
gluon channel. 

4.4 The Spatial Wilson Loop 

The spatial Wilson loop W in 3D is not a static operator, and thus not predicted to 
assume the same value at high temperature and in the 2D reduced model. However i) 
a non zero spatial string tension &2+1 is known to exist |2(J above the deconfinement 
temperature, ii) the pure gauge 2D theory is confining and produces a string tension 
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Figure 5: Physical screening masses Ms in units of g^vT versus gf/T, in (2 + 1)D 
(black points) and 2D (squares). Also shown are the masses obtained with the 
numerical value of h 2 from Eq. ( p6|) (triangles) instead of its asymptotic expression 
(|25|), and with the "naive" reduced action: no Higgs potential, h 2 = h± = 0. 



cr" and iii) in two dimensions we expect confinement and thus a finite string 
tension a 2 to survive when the Higgs field <f> is turned on. 

It is thus interesting to compare these three quantities as a function of the 
temperature of the (2 + 1)D model. For this comparison, we take a 2 o 2 +\ and o?o\ 
from the analysis of W data provided by [00 and by our simulation, whereas a 2 a 2 
is computed analytically at large (3 2 = L /3s. This is obtained by expanding in l//? 2 
the one plaquette SU(3) partition function (Eq. (^) for one single site). The SU(3) 
integral is performed following p0[ ; an d one finds 



aV 2 = -i + -^ + O(/? 2 - 3 ). (53) 
The next terms are easy to derive, but not required for our purpose. This expression 
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Figure 6: The square root of the physical string tension in units of g^\T as a 
function of g\jT in (2+l)D (filled circles) and 2D (squares). The line denotes the 
scaling limit a/2/3 in 2D pure gauge theory (|54 



can be equivalently rewritten 



-5^- = - + — — aT 2 + O 

g\T 3 36 T ; 



(54) 



showing that a\ scales as 2g|T/3, up to scaling violations at finite T of order {aT) 2 = 
1/Z/q. The quantities reported in Fig. 6 versus g\jT are the values of yfa in units 
of g3\/T, which we thus compare to a/2/3. The numerical values of ^/{glT) and 
cr 2+i/ (#3^) were obtained from the simulations for Lo=4. From the 2D point of view, 
we see that introducing the 4> field modifies the picture quite weakly: no sizable g\jT 
dependence around a value close to a/2/3. The behaviour for a 2 +\ is substantially 
different, with a sizable slope in g\jT . The difference observed between a 2 +i and a\ 
may be understood as a consequence of W not being a static operator: the average 
of its non static modes with the (2 + 1)D weight are missing in its calculation with 
the effective action. However, scaling violations may also contribute differently to 
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these two ex's, so that definite conclusions require complementary simulations. Let 
us finally recall that systematic errors can be present in a measurements, so that 
we cannot infer from the data presented in Fig. 6 whether or not their infinite T 
limit is the same, and equal to 2/3. 



5 Conclusion 



In this article we have described a detailed investigation of the validity of dimensional 
reduction for SU(3) gauge theory in 2+1 dimensions at high temperature. We have 
constructed the reduced model for the static variables in the same way as it was 
previously done in || [7], [8] on the reduction of QCD in 3+1 dimensions. This means 
that the integration over the non static modes has been performed perturbatively, 
and the effective couplings have been kept to one loop order for the two point and 
four point functions. Higher order contributions go to zero at high temperature. 
Higher derivative couplings are unimportant at large distances. 

We investigated the validity of this approximation by calculating the correlation 
function between Polyakov loops, and the spatial string tension. For the 2+1 di- 
mensional theory and the corresponding 2 dimensional reduced adjoint Higgs-gauge 
model we were able to obtain very precise numerical results. Therefore we could 
make a detailed comparison, including a scaling analysis which strongly supports 
the assumption that our lattices correspond to a large enough value of the temporal 
extent L . 

We have found that for the correlation between Polyakov loops the dimensional 
reduction works very well down to temperatures T « 1.5T C , even at distances down 
to or below 1/T. And it does even better than in the 3+1 dimensional case. 

The correlation between Polyakov loops is well described by a simple pole in 
momentum space, even at high temperatures. This is in some contrast to 3+1 
dimensional QCD, where at high temperatures the data seem to indicate the presence 
of a cut, corresponding to the exchange of two Debye screened gluons p3l pTl E4L [34 



The difference may be related to the stronger infrared divergences in 2+1 dimensions, 
making the simple perturbatively resummed Debye screening picture invalid. A 
further investigation of this behaviour in 3+1 dimensions, as well as the investigation 
of other correlation functions in the lower dimensional case would certainly be very 
interesting. 

The spatial Wilson loops in 2+1 dimensions do not correspond to a static opera- 
tor. However, for the string tension, which is extracted from Wilson loops of extent 
larger than 1/T, one can hope that the non static corrections are small. That this 
is in fact the case is supported by the actual comparison with the two dimensional 
model. For this operator, the Higgs sector seems to have little influence, and there 
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is a fairly good agreement also with the pure two dimensional gauge theory, which 
is analytically solvable, and where confinement is given by the two dimensional 
Coulomb potential. The differences between the string tension in the three cases 
considered is a few percent. To judge if these differences are real continuum effects, 
one must further study finite size and scaling corrections. 

Our work shows that it may be possible to explain the non perturbative fea- 
tures of 2+1 dimensional QCD in the deconfined phase with a relatively simple two 
dimensional model. 



Acknowledgments The authors thank Christian Legeland for the use of his 
3D SU(3) data and are grateful to Jean-Paul Blaizot for his interest and clarifying 
comments. P.B. was supported by the Alexander von Humboldt Foundation and 
partially by the KBN grant 2 P03B 019 17, and T.R. by a Heisenberg fellowship. K.P. 
was supported by DAAD. We also thank the DFG for support under the contract Ka 
1198/4-1. This project was carried out in part on the new cluster computer ALICE 
at COMPASS, Wuppertal University, during its initial test phase. We thank the 
staff of the ALICE cluster-computer at Compass, University of Wuppertal, for their 
support. 

A Calculation of the (p 2 coupling in the Effective 
2D Action 



We derive the result for hi announced in Eq. (p5|) . According to Eqs.(^g, [U| p0|), 
the sum we need to perform in the large Lq limit is 

^e-^(0) = \y d ^ x = ^ (55) 
3 oov ) 2 ^ dx , L , v I 

no=l 

G(x) = sin(x) cos(x) S'(sin 2 (x)) (56) 

1 r dkl dkl 

[t} = 4^ 7_ 7r 2(1 + t) - cos(fci) - cos(k 2 ) ' ( °' ) 

In the latter integral, we replace the inverse denominator D^ 1 by J °° dy exp(—y D) , 
integrate over k±, ki and use the definition of the modified Bessel function Iq to get 



POO 

S(t)= / dyex V {-2(l+t)y]lZ(y). (5? 
Jo 
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A closed form for this integral is (see for example Eqs. (6.612.4) in |3l[] and (17.3.9) 
in H|) 

-zf/ 2 , 1 



S(t) 



-F(l/2,l/2;l;z), z 



(59) 



which relates the behaviour of S(t) at low t to that of the hypergeometric function 
F(l/2, 1/2; 1; z) at large negative z. One finds (Eq.(15.3.13) of Ref. 



1 oo 
n=0 



r(n + l/2) 

r(n + i) 



log(-z) + 2m(n + 1) - 2m(n + 1/2) 



(60) 



A first approximation Wq to W in (55) is obtained from the standard relationship 
between a sum and an integral: 



La 
2tt 



tt — it/Lo i 

dG(x) + - 

tt/Lo 4 



~dG(x) 


dG(x) 




. dx 




dx 


X = TT--P-. 

L 



L o^f n \ , ldG ( x ) 
it Lq 2 dx 



(61) 
(62) 



where use has been made of the symmetry of under x <-> 7r — x. Up to terms 
which vanish as vt/Lq — > 0, Eq. fl5U| ) gives: 



3 log 2 



+ 1 - log 7T + log L c 



(63) 



For that part of G(x) which is analytic at x = 0, this is the right answer. A 
correction to the constant term in Lq however arises from the logarithmic singularity 

of G{x) ^Euler — Mac-Laurin formula, As x — > 0, we have ^ ~ — ~log^ 

up to analytic terms, whose contribution W l ° 9 to W from Eq. 
using 

L /2 



710 = 1 



Y] log 7^ = lo 

^0 



(Wg) 

T L /2 



log 



(2e) L °/ 2 



can be computed 
(64) 



and found to be 



W l ° 9 = 

2tt 



log(2L o7 r) -L (l + log2) 



(65) 



This we compare to the contribution Wq° 9 to Wo of the same singular part of ^P, 
namely 



Ztt 



log(2L )-L (l + log2) 



(66) 
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so that in the limit L { 



'0 



oo the final result is 



W = W -Wq O9 + W 1 ° 9 

= -2^ 1Og2 - 1 + 1OgL0 ] 



(67) 
(68) 
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